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Executive  Summary 

This  work  involved  the  study  of  instabilities  and  transition  to  turbulence  in  boundary 
layers  on  cones  at  different  angles  of  attacks  at  high  supersonic  and  hypersonic  Mach  num¬ 
bers.  It  utilized  facilities  at  NASA  Langley  Research  Center  which  make  use  of  “quiet” 
design  nozzles  to  minimize  the  level  of  acoustic  disturbances.  The  experiments  involved  two 
models.  The  first  was  a  7°  half-angle  circular  cross-section  cone  designed  by  Cavalieri[l]. 
It  was  used  in  experiments  at  zero  angle  of  attack  [2,  3],  where  the  principle  instability  is 
pairs  of  traveling  oblique  wave  pairs;  and  at  angles  of  attack  where  the  principle  instability 
is  cross-flow  modes.  The  circular  cone  at  an  angle  of  attack  has  a  3-D  boundary  layer  which 
is  similar  to  that  on  elliptic  cone,  which  was  the  second  model  built.  This  consisted  of  a  10° 
half-angle  cone  with  a  4:1  elliptic  cross-section.  The  elliptic  cone  was  designed  to  operate 
in  a  new  Mach  8  tunnel  at  NASA  Langley.  It  mounts  on  a  sting  which  placed  it  inside  the 
worst-case  estimate  of  the  location  of  the  nozzle  quiet  zone.  The  model  is  annealed  and 
hardened  to  withstand  the  900°F  free-stream  temperature,  and  maintain  a  polished  finish. 
A  3-D  traversing  mechanism  capable  of  following  the  tapered  elliptic  surface,  and  operating 
at  high  temperatures  in  the  hypersonic  flow  was  also  designed  [4].  Unfortunately,  the  elliptic 
■£one  was  never  used  because  of  technical  problems  with  the  Mach  8  tunnel  which  prevented 
%  from  ever  operating.  However,  the  elliptic  cone  was  also  designed  to  operate  in  the  Mach 
3.5  tunnel  used  for  the  circular  cone  experiments.  A  pending  AFOSR  proposal  describes  an 
experimental  program  for  its  use  at  the  lower  Mach  number. 

The  research  program  focussed  on  the  introduction  of  controlled  disturbances  to  excite 
the  most  unstable  linear  instability  modes,  and  on  quantitative  time-resolved  measurements 
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of  the  instabilities  for  direct  comparison  to  numerical  simulations.  Both  of  these  were  suc¬ 
cessfully  completed.  The  controlled  disturbances  were  based  on  phased  plasma  actuators. 
During  the  course  of  the  research,  these  were  perfected  to  operate  for  virtually  infinite  time 
without  degradation.  The  time  resolved  measurements  were  made  using  hot-wires.  For  these 
we  developed  an  in-situ  calibration  method  to  obtain  the  mean  quantities,  (pU)  and  T0,  and 
the  fluctuating  quantities  ( pU )'2  and  Tq2.  This  was  tested  with  the  sharp-tipped  circular  cone 
at  zero  angle  of  attack  in  order  to  quantify  the  error  in  in-situ  hot-wire  calibration.  These 
measurements  also  led  to  important  parameters  in  the  probe  geometry  that  were  needed  to 
minimize  its  upstream  influence  in  the  subsonic  portion  of  the  boundary  layer. 


Contents 


1  Background  1 

2  Experimental  Approach  2 

2.1  Cone  Models  .  2 

2.2  Plasma  Actuator  Array .  3 

2.3  Results:  Circular  Cone  Actuation .  6 

2.4  In-situ  Hot-Wire  Calibration .  8 

3  Personnel  and  Technology  Transfer  16 

4  Summary  16 


1  Background 

In  recent  years  there  has  been  a  renewed  interest  in  clarifying  the  mechanisms  of  stability  and 
transition  to  turbulence  of  high  speed  compressible  boundary  layers.  Improved  predictions 
of  transition  Reynolds  numbers  are  needed  in  aeronautics  applications  such  as  high-speed 
vehicles  or  Earth  re-entry  systems.  This  motivation  reflects  the  current  evolution  of  the 
US  Air  Force  from  its  static  Cold  War  posture  to  one  that  is  capable  of  projecting  mili¬ 
tary  power  anywhere  in  the  world  at  short  notice,  i.e.  Global  Reach,  Global  Power.  The 
ability  to  formulate  predictions  of  transition  Reynolds  numbers  for  high-speed  flows  on  3-D 
lifting  bodies  requires  a  fundamental  understanding  of  the  roles  of  instability  mechanisms 
in  transition  behavior.  In  spite  of  its  importance,  there  is  very  little  accurate  data  with 
which  to  develop  empirical  criteria  or  to  validate  numerical  simulations  that  might  be  used 
in  design.  To  obtain  this  will  require  a  concerted  effort  which  involves  a  partnership  between 
experimental  and  computational  investigations.  The  challenge  for  experiments  is  to  control 
unsteady  initial  conditions  so  that  in  direct  comparisons  to  simulations,  initial  amplitudes 
and  wave  numbers  can  be  specified.  The  outgrowth  of  such  control  is  the  potential  to  modify 
the  transition  front  and  therefore  define  initial  conditions  which  can  enhance  performance 
by  reducing  drag  and  heat  loading. 

In  all  the  past  experiments,  except  those  of  Kosinov  et  al.[5]  transition  to  turbulence 
developed  from  background  disturbances.  The  problem  with  that  approach  is  that  the  char¬ 
acter  of  the  free-stream  disturbances  (both  vortical  and  acoustic)  are  difficult  to  completely 
document,  and  are  largely  facility  dependent.  In  supersonic  tunnels  with  free-stream  Mach 
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numbers  greater  than  2,  disturbances  originate  from  side-wall  boundary  layers  which  are 
unstable  and  become  turbulent.  These  are  felt  in  the  center  of  the  nozzle  as  acoustic  (Mach) 
waves.  The  “quiet”  design  supersonic  nozzles  [7]  overcome  this  through  a  combination  of 
wall  suction  at  the  throat,  and  well  designed  nozzle  expansion  sections  which  have  curvature 
which  maximizes  the  transition  Reynolds  number  for  a  Goertler  instability. 

Even  if  the  background  disturbances  are  well  documented,  a  second  problem  is  that  the 
manner  by  which  free-stream  disturbances  enter  and  excite  boundary  layer  modes  (recep¬ 
tivity)  is  dependent  on  the  model  geometry  such  as  local  curvature,  surface  junctions,  or 
imperfections  (micro  roughness). 

Kosinov  et  al.[5]  aimed  at  overcoming  these  problems  by  introducing  periodic  distur¬ 
bances  at  a  point  on  the  surface  of  a  flat  plate  at  M  =  2.0.  The  disturbances  were  produced 
by  a  glow-discharge  between  two  electrodes  in  a  small  cavity  below  the  plate  surface.  The 
forcing  at  a  point  produced  a  full  spectrum  of  modes  with  different  wave  angles.  Ultimately 
the  most  amplified  mode  was  determined  by  the  flow  condition,  and  was  in  this  regard, 
uncontrolled. 


2  Experimental  Approach 

In  this  work,  we  have  tried  to  develop  an  experimental  analog  to  DNS  calculations  whereby 
we  could  specify  the  amplitudes,  frequency,  and  wave  angles  of  disturbances.  The  approach 
follows  that  used  in  numerous  incompressible  flows  [6,  8]  with  the  exception  of  the  actuator. 
In  order  to  minimize  the  required  levels  of  controlled  disturbances,  the  experiment  was 
performed  in  the  0.5  m.  Mach  3.5  “quiet”  tunnel  at  NASA  Langley  Research  Center. 


2.1  Cone  Models 

The  first  model  was  a  <j)c  =  7°  circular  cone  that  is  shown  in  Figure  1.  The  cone  is  35.56  cm. 
in  length.  Its  length  was  chosen  to  be  compatible  with  the  Mach  3.5  and  61  “quiet”  tunnels 
at  NASA  Langley.  At  the  typical  stagnation  pressures  used  in  the  past[3],  the  first  30.48  cm 
(86%)  of  the  model  is  inside  the  “quiet  zone.” 

The  circular  cone  is  hollow  with  thin  walls.  It  is  instrumented  with  streamwise  arrays  of 
thermocouples  and  pressure  ports.  Wiring  for  the  actuator  array  also  run  through  the  cone 
to  an  actuator  hub  which  is  14.3  cm.  from  the  cone  tip.  The  cone  tip  is  removable  to  allow 
different  tip  bluntness.  The  “sharp”  tip  has  a  nose  radius  of  0.038  mm. 

The  circular  cone  is  supported  by  a  sting  which  also  holds  a  motorized  3-D  traversing 
mechanism.  The  traversing  mechanism  holds  removable  wedges  in  which  miniature  pitot  or 
hot-wire  probes  are  held.  Further  details  of  the  circular  cone  and  traversing  mechanism  are 
given  by  Cavalieri[l]. 

The  other  model  is  a  4:1  elliptic  cross-section  cone.  A  photograph  is  shown  in  Figure 
2.  The  angles  of  the  major  and  minor  axes  are  <f)c  =  10°  and  2.65°,  respectively.  This  cone 
was  originally  designed  for  a  Langley  Mach  8  “quiet”  tunnel  which  was  never  completed.  It 
however  is  the  same  length  (35.56  cm)  as  the  circular  cone,  and  will  fit  within  the  quiet  zone 
in  the  Langley  Mach  3.5  tunnel. 

1This  model  was  used  once  in  this  tunnel  before  it  was  dismantled  at  Langley.  Efforts  are  being  made  to 
reopen  this  tunnel. 
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Figure  1:  Photograph  of  circular  cone  and  traversing  Figure  2:  Photograph  of  the  4:1  el- 
mechanism  shown  outside  of  the  0.5  meter  Mach  3.5  liptic  cone  and  support  sting, 
“quiet”  tunnel  at  NASA  Langley. 


Figure  3:  Illustration  of  traversing  mechanism  design  for  the  elliptic  cone. 


The  elliptic  cone  and  sting  are  fabricated  from  416  stainless  steel  and  annealed  and 
hardened  to  withstand  up  to  1360°R  (the  design  stagnation  temperature  of  the  previously 
proposed  Mach  8  tunnel).  A  design  for  a  traversing  mechanism  has  been  developed  by 
Siddiqi[4j.  It  is  shown  in  the  illustration  in  Figure  3. 

2.2  Plasma  Actuator  Array 

Our  objective  was  to  have  actuators  that  could  introduce  periodic  disturbances  that  would 
excite  instability  modes  consisting  of  pairs  of  helical  waves  at  the  same  frequency  and  equal- 
opposite  wave  angles. 

Our  approach  was  based  on  using  an  azimuthal  array  of  disturbance  actuators  at  one 
streamwise  location  [8].  Critical  aspects  of  the  design  were  the  azimuthal  spacing,  s,  of 
the  array  elements,  and  the  streamwise  location  of  the  actuator  on  the  cone.  The  ideal 
streamwise  position  of  the  actuator  is  near  the  location  of  the  lower  neutral  branch  of  the 
linear  stability  growth  curve.  The  choice  of  s  is  based  on  the  range  of  frequencies  and  wave 
angles  of  the  most  amplified  modes. 

Linear  stability  calculations  indicated  that  the  most  amplified  dimensionless  frequencies 
were  in  the  range  from  10  <  F  x  106  <  30,  and  had  wave  angles,  ip,  of  approximately  60°. 
In  order  to  maximize  the  extent  of  the  quiet  zone,  the  experiment  was  designed  to  operate 
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Table  1:  Excitation  Conditions  based  on  an  azimuthal  spacing  of  15° 


f  (kHz) 

F  x  106 

<M°) 

71  —  27TF actuator  / ^0 

50 

11.9 

66 

~  12 

60 

14.2 

62 

12 

70 

16.7 

58 

12 

Figure  4:  Cut-away  illustration  of  plasma  actuator  Figure  5:  Film  clip  of  new  actuator 

hub  used  to  excite  T-S  modes  on  circular  cylinder.  operating  in  a  low  pressure  cham¬ 

ber  with  180°  phase  shift  between 
anodes. 

at  the  lowest  possible  stagnation  pressures.  We  chose  Pq  to  be  172.35  kPa  (25  psi.)  At  this 
value,  and  with  To  =  311K,  the  actuator  location  of  x  =  14.303  cm.  placed  it  just  upstream 
of  the  lower  neutral  branch  for  the  band  of  most  amplified  frequencies. 

We  chose  an  azimuthal  spacing  of  15°  for  the  individual  array  elements.  Based  on  this, 
Table  1  gives  the  largest  wave  angles  that  were  possible  for  the  range  of  frequencies  used  in 
the  experiment. 

The  actuator  is  based  on  an  a.c.  glow  discharge.  A  schematic  cut-away  drawing  of  the 
actuator  hub  is  shown  in  Figure  4.  The  actuator  is  made  up  of  azimuthally  spaced  anodes 
that  are  supplied  with  a  high-voltage  a.c.  signal.  The  frequency  of  the  a.c.  matches  the 
instability  excitation  frequency.  The  relative  phase  of  the  a.c.  between  adjacent  anodes 
controls  the  instability  wave  angle.  The  metal  tip  of  the  cone  is  grounded  and  therefore 
forms  a  common  cathode.  When  the  voltage  is  high  enough,  ionized  air  (plasma)  forms  in 
the  cavity.  This  results  in  a  body  force  on  the  non-ionized  air,  which  induces  an  unsteady 
jet  of  fluid  that  exits  through  the  hole  in  the  cone  surface. 

Figure  5  shows  two  frames  of  a  film  clip  of  the  actuator  operating  in  a  low  pressure 
chamber,  which  simulates  the  static  pressure  in  the  Mach  3.5  flow  over  the  cone.  In  this 
case,  the  individual  anodes  were  operated  with  a  180°  phase  shift.  The  plasma  in  the 
photographs  appears  blue. 
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Figure  6:  Photographs  of  two  plasma  actuators  which  were  designed  to  excite  different  wavelengths 
of  cross-flow  modes  on  subsonic  (left)  and  supersonic  (right)  swept  wings. 


A  different  plasma  actuator  has  been  designed  for  principally  exciting  stationary  cross- 
flow  modes.  This  would  be  used  with  the  cylinder  at  an  angle  of  attack,  or  with  the  elliptic 
cone. 

The  concept  is  based  on  the  work  of  Saric  et  al.  [9],  in  which  the  transition  Reynolds 
number  of  3-D  boundary  layers  subject  to  a  cross-flow  instability,  could  be  substantially 
increased  by  surface  bumps  which  excite  less  amplified  spanwise  wave  lengths.  In  association 
with  Saric  in  joint  DARPA  funded  research,  we  developed  plasma  actuators  which  could 
produce  an  effect  equivalent  to  the  surface  bumps. 

The  actuators  are  fabricated  on  0.1  mm  thick  F-4  grade  fiberglass  re-enforced  epoxy  strips 
which  are  copper  clad  on  both  sides  (0.025  mm  thick  copper  clad  Kapton  is  also  available). 
The  copper  is  photo-etched  to  leave  a  pattern  of  conducting  lines  that  serve  as  the  anodes 
and  cathodes.  The  surfaces  are  subsequently  doped  to  leave  a  smooth  uniform  surface  that 
is  sanded  smooth,  with  the  electrodes  flush  to  the  surface. 

The  plasma  forms  at  the  edges  of  the  anodes  where  they  overlap  the  cathode.  The 
body  force  is  such  that  it  induces  a  flow  away  from  these  edges.  This  results  in  a  thinning 
of  the  boundary  layer  over  the  anodes,  and  a  thickening  of  the  boundary  layer  between 
the  anodes.  The  wavelength  of  the  alternate  thickening  and  thinning  is  specified  by  the 
distance  between  anodes.  The  two  actuators  were  designed  for  two  swept  wing  experiments 
at  ASU;  for  subsonic  (left)  and  supersonic  (right)  wings.  The  advantage  of  the  “plasma 
bumps”  over  fixed  surface  bumps  is  that  the  spacing  can  be  modified  to  suit  changes  in  the 
flight  conditions.  In  addition,  the  plasma  actuator  can  be  operated  to  produce  traveling 
disturbances  which  in  the  3-D  boundary  layer,  would  excite  traveling  cross-flow  modes. 
These  are  more  amplified  than  the  stationary  modes,  and  would  be  preferred  if  a  more  rapid 
transition  is  needed.  Finally,  if  conditions  warrant,  the  actuator  can  be  turned  off,  which  is 
not  possible  with  fixed  surface  bumps. 
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Figure  7:  Spectra  at  different  streamwise  locations  in  the  boundary  layer,  at  a  fixed  height  of 
y  =  0.356  mm.,  with  excitation  at  50kHz  (F  =  11.9,  ip  =  ±66°). 


2.3  Results:  Circular  Cone  Actuation 

All  of  the  controlled  instability  results  were  for  the  circular  cone  at  zero  angle  of  attack.  These 
used  the  actuator  design  which  was  illustrated  in  Figure  4.  The  excitation  was  designed  to 
excite  pairs  of  helical  modes  with  azimuthal  wave  angles  of  ±66°.  The  excitation  frequency 
was  50kHz,  which  corresponds  to  the  conditions  in  the  top  row  of  Table  1. 

An  example  of  the  results  with  excitation  are  shown  in  Figure  7.  The  figure  shows  spectra 
at  a  fixed  height  in  the  boundary  layer,  and  different  downstream  distances  (15.24  <  x  < 
27.94  cm).  The  physical  frequency  is  nondimensionalized  as  F  —  2-nfv/lf1. 

The  peaks  at  F  =  11.9  and  23.8  correspond  to  the  excitation  frequency  (fe)  and  its 
harmonic  (2/e),  respectively.  The  most  interesting  feature  of  these  spectra  is  the  growth  of  a 
subharmonic  of  the  excitation  frequency  (fe/ 2)  at  the  further  downstream  positions.  Along 
with  the  subharmonic  peak,  a  peak  at  3/e/2  is  also  observed.  This  likely  evolves  from  a 
summing  interaction  between  the  fundamental  and  subharmonic  modes. 

In  order  to  better  understand  the  origin  of  the  subharmonic,  the  effect  of  changing  the 
initial  50kHz  amplitude  was  investigated.  This  is  shown  in  the  top  part  of  Figure  8.  This 
was  documented  at  x  =  23.5  cm.,  y  =  0.356  mm.  and  6  —  0°.  Here  the  initial  amplitude 
above  a  minimum  to  produce  a  plasma,  was  varied  by  a  factor  of  approximately  two.  The 
abscissa  corresponds  to  the  measured  amplitude  at  the  fundamental  frequency  normalized 
by  its  minimum  value.  The  accuracy  of  the  amplitude  is  0.05  percent  of  the  maximum  value. 

The  ordinate  in  the  plot  is  the  normalized  amplitude  at  the  subharmonic  frequency  on  a 
log  axis.  If  the  subharmonic  mode  results  from  a  secondary  interaction  with  the  fundamental, 
we  expect  it  to  have  an  exponential  response  to  the  fundamental  mode  amplitude.  The  results 
seem  to  support  this,  since  the  points  fall  on  a  straight  line  on  a  log  axis. 
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Figure  8:  Response  of  the  subharmonic  mode  to  different  fundamental  mode  amplitudes 
(top), and  subharmonic  amplitude  azimuthal  variation  (50kHz  excitation  measured  at  x  = 
23.5  cm.,  y  =  0.356  mm.). 
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Measurements  were  taken  at  different  azimuthal  positions  to  determine  the  azimuthal 
wave  length  of  the  subharmonic  mode.  The  result  is  shown  in  the  lower  plot  in  Figure  8. 
This  corresponds  to  the  same  spatial  position  as  before.  The  azimuthal  angle  ranges  from 
±25°  with  respect  to  the  0°  reference  position  used  in  the  previous  figures.  The  accuracy 
in  the  azimuthal  position  is  0.56°.  The  subharmonic  amplitude  has  been  normalized  by  the 
value  at  0  =  0°. 

The  azimuthal  amplitude  variation  for  helical  mode  pairs  should  consist  of  periodic  max¬ 
ima  and  minima.  For  the  n  =  12  fundamental  mode  that  was  excited,  amplitude  minima 
should  occur  approximately  30°  apart.  A  secondary  instability  analysis  in  the  boundary 
layer  on  a  flat  plate  at  Mach  1.6  by  Ng  and  Zang7  which  specified  a  primary  instability  of 
pairs  of  oblique  waves  at  a  fundamental  frequency,  found  the  most  amplified  subharmonic 
consisted  of  pairs  of  oblique  waves  of  comparable  wave  angles.  This  appears  to  be  consistent 
with  the  spacing  of  the  minima  for  the  subharmonic  mode  in  our  experiment. 

The  results  in  the  previous  two  figures  were  all  based  on  time  series  obtained  from  a 
hot-wire  which  was  uncalibrated,  but  operated  at  a  fixed  high  overheat  so  that  it  should 
be  primarily  sensitive  to  mass-flux  fluctuations.  This  is  a  limitation  of  these  measurements 
which  we  address  with  a  new  in-situ  calibration  approach  that  was  developed  as  part  of  this 
grant. 


2.4  In-situ  Hot-Wire  Calibration 

To  our  knowledge,  there  has  never  been  a  successful  attempt  to  measure  amplification  rates 
with  calibrated  how-wires  in  supersonic  boundary  layers.  Typically,  un-calibrated  wires 
are  used,  with  the  assumption  that  when  maintained  at  a  high  enough  overheat,  the  wire 
becomes  mostly  sensitive  to  mass-flux  fluctuations  and  the  amplification  rate  of  linear  modes 
can  be  deduced  by  comparing  voltage  fluctuation  magnitudes  at  different  spatial  locations. 
With  un-calibrated  sensors,  this  requires  that  at  every  measurement  station,  the  wire  is 
placed  at  a  point  within  the  boundary  layer  where  the  mean  mass-flux  is  constant  to  ensure 
that  the  dynamic  sensitivity  remains  fixed.  The  difficulties  associated  from  trying  to  infer 
amplification  rates  with  un-calibrated  sensors  becomes  more  severe  when  trying  to  measure 
nonlinear  disturbances,  such  as  the  development  of  sub-harmonic  resonance  that  we  observed 
[3]  in  Figure  7. 

The  appearance  of  these  resonant  modes,  resulting  from  non-linear  interactions  with  the 
linear  modes,  stresses  the  need  for  calibrated  time-resolved  measurements.  Our  efforts  in  the 
last  year  of  the  grant  had  therefore  been  directed  towards  developing  an  “in-situ”  approach 
to  hot-wire  calibration,  which  is  described  in  the  following. 

The  standard  approach  to  calibration  involves  placing  the  hot-wire  in  the  free-stream  at 
a  fixed  Mach  number  while  changing  stagnation  pressure  ( Pa )  and  temperature  (T0)  (hence 
mass  flux)  at  multiple  overheats.  However,  this  may  not  be  feasible  in  all  facilities  with 
limited  control  over  P0  and  T0.  We  should  also  stress  that  many  calibration  data  pairs  are 
needed  to  achieve  an  accurate  calibration  which  is  suitable  for  measurement  of  fluctuation 
amplitudes. 

Consequently,  an  approach  has  been  developed  that  simultaneously  calibrates  a  hot-wire 
while  performing  wall-normal  surveys  of  the  fluctuating  flow-field.  This  approach  is  termed 
“in-situ”  calibration,  and  makes  use  of  the  variation  of  Pa  and  Ta  within  the  boundary 
layer.  This  is  done  by  performing  total-pressure  surveys  with  a  miniature  probe  to  obtain 
Mach  number  profiles.  Subsequently,  a  hot-wire  is  traversed  through  the  layer  at  multiple 
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Pitot  Probe 


Variable  Overheat 
Hot-wire 


Figure  9:  Schematic  of  process  for  in-situ  hot-wire  calibration. 


Figure  10:  Schematic  for  determining  local  Mach  number  profile  in  the  boundary  layer. 


overheats.  At  the  lowest  overheat,  the  wire  measures  the  stagnation  temperature,  which 
together  with  the  previously  established  Mach  profile  provides  the  mean  mass-flux  and  Ta 
required  for  the  calibration.  This  is  illustrated  in  the  schematic  in  Figure  9. 

This  calibration  technique  was  carried  out  in  the  boundary  layer  over  the  circular  cone 
at  zero  angle  of  attack  at  Mach  3.5.  This  case  was  selected  because  we  could  easily  solve  the 
mean  flow  similarity  solution  for  this  conic  flow  and  determine  the  error  in  our  measurements. 
It  begins  with  the  measurement  of  the  Mach  number  profile  at  an  rr-station  in  the  boundary 
layer.  For  this,  the  static  pressure  is  measured  using  a  wall  pressure  port  in  the  model,  and 
the  total  pressure  is  measured  with  a  miniature  pitot  probe  that  is  traversed  through  the 
layer.  This  is  illustrated  in  Figure  10. 

The  Mach  number  profile  is  found  from  a  recursive  solution  to  the  identity 


p  _  ^02  P*  p 

^  5  -A 

F(M)  HM) 
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This  step  represents  the  upper-left  box  in  Figure  9. 

This  profile  is  repeated  at  the  same  heights  above  the  wall  with  a  hot-wire  which  is 
operated  at  different  (10  or  more)  over-heat  values.  The  stagnation  temperature,  T0,  is 
found  from  the  hot-wire  at  zero  overheat  (cold-wire).  Based  on  the  resistive- temperature 
properties  of  the  wire  material,  the  wire  temperature  is 


(■♦*£>)• 


Tw  —  Tref  + 

Then  the  stagnation  temperature  is 

To  =  Tw/r) 

where  77  is  the  recovery  factor. 

Knowing  To,  the  static  temperature  is  found  from 

m  T0 


1  +  ifM2' 


Finally,  the  mass  flux  profile  is  found  from 

(pUj  =  p°iM  ^ 


yfT  V  R‘ 


(1) 

(2) 

(3) 

(4) 


Implicit  in  the  in-situ  calibration  approach  is  the  need  to  accurately  measure  the  mean 
quantities  through  the  boundary  layer.  A  somewhat  unexpected  issue  that  emerged  in  this 
was  the  effect  of  probe  interference.  An  example  is  shown  in  Figure  11.  Here  we  observed 
two  deviations  in  the  measured  values  from  the  theoretical  distribution:  an  over-shoot  at 
the  edge  of  the  boundary  layer,  and  an  under-shoot  in  the  lower  part  of  the  layer. 

These  two  types  of  experimental  deviations  had  previously  been  observed  by  Kendall[10] 
and  Morkovin  and  Bradfield[l  1] .  In  both  cases  it  was  associated  with  probe  interference.  At 
the  outer  edge  of  the  boundary  layer,  Kendall  could  eliminate  the  overshoot  when  the  pitot 
probe  diameter  was  no  larger  than  12  percent  of  the  boundary  layer  thickness.  Morkovin 
and  Bradfield  associated  the  deviation  in  the  lower  part  of  the  boundary  layer  with  an 
upstream  rise  in  the  static  pressure  produced  by  the  sensor  and  communicated  upstream 
in  the  subsonic  portion  of  the  boundary  layer.  They  recommended  an  upstream  influence 
distance,  s,  where  the  pressure  disturbance  diminishes  by  e-1  given  as 


s  = 


1.3  Tw 
VM2  -  1  Te 


M2  +  2  (Tw\07 

2s/M2  -  1  VtJ 


(5) 


For  our  experimental  conditions,  s  ~  105. 

The  recommendations  from  these  early  investigations  was  used  to  obtain  the  results 
shown  in  Figure  12.  In  this  case,  the  pitot  tube  exterior  dimensions  were  0.203  mm.  wide 
by  0.102  mm.  high  (0.008  in.  wide  by  0.004  in.  high).  For  a  nominal  boundary  layer  thick¬ 
ness,  the  probe  height  was  approximately  0.085,  which  satisfied  Kendall’s[10]  guideline.  The 
probe  extended  1.27  cm.  (0.5  in.)  upstream  of  the  probe  support  wedge.  This  corresponds 
to  approximately  105  in  accordance  with  the  Morkovin  and  Bradfield[ll]  guidelines.  A  pho¬ 
tograph  of  the  final  arrangement  is  shown  in  Figure  13.  The  tubulation  used  is  designated 
as  “e”. 
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Figure  11:  Mach  number  profile  at  one  x- 
location  in  the  boundary  layer  over  the  cir¬ 
cular  cone.  M  =  3.5,  a  =  0.  Solid  curve  is 
theoretical  profile. 


Figure  12:  Mach  number  profile  at  one  re¬ 
location  in  the  boundary  layer  over  the  cir¬ 
cular  cone  measured  with  minimum  probe  in¬ 
terference.  M  =  3.5,  a  =  0.  Solid  curve  is 
theoretical  profile. 


a:  0.013  x  0.007  in  d:  0.008  x  0.004  in 

b:  0.013  x  0.007  in  (squashed)  e:  0.008  x  0.004  in  (squashed) 

c:  0.010  x  0.005  in 


Figure  13:  Photographs  of  probe  support  wedge  with  pitot  probe  at  tip,  and  cross-sections  of 
tubulations  used  to  measure  Mach  number  profiles.  Results  in  Figure  12  used  tubing  “e” . 
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Figure  14:  Comparison  of  the  cold- wire  tem¬ 
perature  through  the  boundary  layer  against 
theoretical  To  profile  (top).  Variation  in 
recovery  factor  through  the  boundary  layer 
(bottom) . 


Mach 


Figure  15:  Recovery  factor  based  on  all  mea¬ 
sured  profiles  over  the  circular  cone.  Curves 
represent  free  molecular  limit  and  Lowell  con¬ 
tinuum  model. 


Figure  12  demonstrates  that  when  following  the  guidelines,  it  was  possible  to  accurately 
measure  the  mean  Mach  number  profiles.  Of  the  two  criteria,  the  Morkovin  and  Bradfield[ll] 
is  the  most  difficult  to  satisfy  since  it  is  a  function  of  local  flow  properties  ( Rx ,  M,TW/Te). 
Therefore  information  about  the  mean  flow  needs  to  be  known  to  some  extent  in  order  to 
properly  design  probes.  This  also  applies  to  the  hot-wire  probes.  We  have  experimentally 
determined  that  when  our  hot-wire  probes  are  designed  by  the  same  criteria  as  the  pito  probe, 
the  upstream  effect,  measured  as  a  local  change  in  the  wall  static  pressure,  is  comparable  to 
the  pitot  probe.  Therefore  we  expect  similar  accuracy  with  the  properly  designed  hot-wire 
probes. 

The  hot-wire,  operating  at  zero  overheat  (cold  wire)  was  used  to  measure  the  T0.  A 
question  that  arises  at  this  point  is  the  value  of  the  recovery  factor,  77.  This  was  investigated 
by  comparing  the  wall-normal  profile  of  Tw  in  the  boundary  layer  against  the  theoretical 
distribution  of  To.  This  is  shown  in  Figure  14.  The  value  of  17  which  most  often  appears 
in  the  literature  is  0.965.  As  these  results  indicate,  the  recovery  factor  is  not  constant,  but 
varies  through  the  layer.  This  was  compiled  for  all  of  the  x  locations  in  Figure  15.  This 
indicates  that  77  varies  from  approximately  1.02  to  0.965.  We  have  shown  that  assuming  a 
constant  value  for  77  in  this  range  accounts  for  at  most  a  3%  error  in  the  mass  flux. 

Knowing  T0  and  M,  Eqs[3  and  4]  are  used  to  determine  (pU).  A  sample  of  the  measure¬ 
ments  obtained  from  the  pitot  probe  is  shown  as  the  square  symbols  in  Figure  16.  The  solid 
curve  is  the  theoretical  profile. 

_  The  approach  to  calibrate  the  hot-wire  is  to  decompose  the  anemometer  mean  voltages, 

V0,  into  (pU)  and  To  through  calibration  for  many  (typically  10)  over-heats.  The  hot-wire  is 
operated  in  a  constant  current  mode  and  frequency  compensated. 

Analysis  of  the  steady  heat  transfer  balance  for  the  hot-wire  operating  in  a  constant 
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Figure  16:  Profile  of  ( pU )  through  the  boundary  layer  for  the  circular  cone  measured  using  the 
pitot  probe  and  cold-wire  (square  symbols),  and  measured  using  the  calibrated  hot-wire  (triangle 
symbols)  Solid  curve  is  the  theoretical  profile. 


current  mode  leads  to  the  following  equation 

cj?Z+c^-c3%1+‘+cMpu)',TZ‘-°'L+c!,(purTZ‘-°,‘-c6(pu)''%1+'‘-‘"-c,Vo  =  o 

(e) 

where  n  and  a  are  known  constants,  and  C\  •  •  •  C7  are  coefficients  to  be  determined  from 
calibration  at  different  over-heats. 

The  problem  is  formulated  to  minimize  the  least-squared  error.  This  leads  to  the  following 
general  matrix,  where  °lt  * ' '  a7i  represent  the  terms  multiplied  by  the  respective  coefficients 
for  the  zth  [(pU),TQ]  data  pair,  and  the  jth  overheat. 


mo* 
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-g 
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With  the  coefficients  determined  from  the  calibration,  the  time  averaged  voltages  at  each 
spatial  point  from  the  hot-wire  can  be  decomposed  into  ( pU )  and  To  by  solving  the  following 
matrix. 


Ej mlj  Zj mijUiij  ••• 

E j  m2jmij  Ej  rn\j 
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Figure  17:  Calibration  quantities  output  using  Eq[6]  versus  input  for  the  hot-wire. 

where  m\j  —  ^ij^o  1  C% ji  C*  C*5jj  ni/ij  and  m§j  • 

Again,  j  represents  the  over-heat  values  used  in  calibration. 

An  example  of  the  calibration  is  shown  in  Figure  17.  This  shows  the  calibration  data 
versus  the  values  obtained  from  the  calibration  relation  (Eq[6]).  The  input-output  relation 

falls  on  a  straight  line  indicating  perfect  agreement.  _ 

A  further  indication  of  the  calibration  comes  from  comparing  the  ( pU )  variation  across 
the  boundary  layer  with  the  theoretical  profile.  This  is  shown  as  the  triangle  symbols  in 
Figure  16.  The  error  in  the  hot-wire  distribution  is  identical  with  that  of  the  pitot  probe. 
Since  one  was  used  to  calibrate  the  other,  this  indicates  a  perfect  calibration  result.  If  the 
theoretical  distribution  had  been  used  to  supply  the  calibration  data  pairs,  the  comparison  to 
theory  could  have  been  perfect.  However,  since  the  theoretical  distribution  may  not  always 
be  known,  we  chose  to  use  calibration  data  from  the  pitot  probe  as  a  complete  test  of  the 
in-situ  calibration  approach. 

Ultimately,  as  illustrated  in  Figure  9,  the  calibrated  hot-wire  is  intended  for  time-resolved 
measurements  such  as  the  variance  of  the  fluctuating  flow  quantities,  ( pU )'2  and  Tq2  in  terms 
of  the  variance  of  the  anemometer  voltage,  V'2. 

In  order  to  relate  these  statistics  to  the  anemometer  voltage  fluctuations,  Eq[6]  is  rewrit¬ 
ten  in  a  more  general  form  as 

Ciiy  v. + a/r;  -  c3T,]  'a  +  cMpurr;-^ + cr>(pijy‘T;;  ~  -  c6(pt/)”r01+“-“  -  c7v„  =  o 

(9) 

For  this,  we  substitute  for  T0,  (pU)  and  VQ  in  terms  of  the  following  decompositions 

T0  =  T0  +  Ti  (10) 

(pU)  =  Q pU )  +  (pU)'  (11) 

V0  =  Vo  +  K'  (12) 
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where  the  primed  terms  represent  the  mean-removed  fluctuations  of  each  of  the  quantities. 

After  substitution,  the  equation  is  expanded  and  linearized.  Within  this  equation,  the 
mean  equation  (Eq[6])  is  identically  zero.  This  leaves  the  following 

qX+q2(pU)'  +q3T^  (13) 

where  _  _ 

<7i  =  Ctf?  +  C4(pU)nT0a~an  -  1  (14) 

q2  =  C,nV0(pU)n~lT^~an  +  C5n(pU)n-lT?~an  -  C6n(pU)n~1%l+a~an  (15) 


q3  =  CiaT0a"Vo  +  C2aT0a  1  -  C3(l  +  a)T?  +  C4(a  -  an)V0(pU)nT{ ?  an  V 

C5(a  -  an)(pU)n7TT~an~l  -  C6(l  +  a  -  an)(pU)n%a~an .  (16) 

Taking  the  variance  of  the  terms  in  Eq[13]  and  rearranging  gives  the  following 

T7=(-l  (pUP+(-)  T?  +  2^Wj%.  (17) 

\qij  \qij  Qiqi 

As  a  general  form,  Eq[17]  can  be  written  as 

k\Di  +  k2D2  +  k3D3  —  k(,  —  0  (18) 

*>=te)2 

which  are  functions  of  the  sensor  overheat,  and 

d1=wp  D2=wrn  d3  =  {£)2t? 

which  are  the  unknown  variances  of  the  flow  quantity  fluctuations. 

The  unknown  quantities  ( Ds )  are  found  by  minimizing  the  least-square-error  defined  as 


a2  =  £  [fcyA  +  k2jD2  +  kZjDz  -  k0jf  (19) 

3= 1 

where  as  before,  j  represents  the  sensor  overheat  values. 

The  unknowns  are  then  found  by  solving  the  following  matrix 


E;  (kfj)  Ej  (k\jk2j)  Ej  {kijhj)  1 

(puy* 

'  -  E j  mojrnij  ' 

Ylj  (&1  j^2j)  Ylj  (^2 j)  Sj  (^2j^3j) 

7^ '2 

i0 

= 

-  E j  m0jm2j 

_  Sj  {kljkzj)  Ylj  (^2j^3j)  Ylj  (^3 j) 

(puyri  J 

.  -  E  j  m0jm3j 
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3  Personnel  and  Technology  Transfer 

This  research  was  performed  by  Professor  Thomas  Corke  along  with  two  graduate  assistants. 
It  benefited  from  the  interaction  of  Mr.  Steve  Wilkinson  at  NASA  Langley  Research  Cen¬ 
ter  who  coordinated  the  Mach  3.5  tunnel  experiments,  as  well  as  made  a  tireless  effort  to 
complete  the  Mach  8  tunnel. 

This  work  will  lead  to  a  Ph.D.  thesis  by  Mr.  Eric  Matlis.  He  is  expected  to  defend 
his  thesis  in  the  December,  2002.  The  other  graduate  student  working  under  this  grant 
developed  the  final  design  of  the  actuator.  She  completed  a  M.S.  thesis[12]  in  August,  2001, 
and  is  presently  pursuing  a  Ph.D.  at  Notre  Dame. 

In  addition  to  the  theses,  one  paper  has  been  published[3],  and  two  are  in  progress  related 
to  the  plasma  actuator,  and  in-situ  hot-wire  calibration  method. 

There  have  been  numerous  spinn-offs  of  the  plasma  actuator  technology  that  was  de¬ 
veloped  under  this  grant  and  a  DARPA  grant.  This  includes  separation  control  on  heli¬ 
copter  rotors  (DARPA/ARO)  and  on  LPT  blades  in  turbo-jet  engines  (NASA  Glenn),  con¬ 
trol  of  jet  flows  to  reduce  sound  (NASA  Langley),  lift  enhancement  (AFIT)[13],  control  of 
cross-flow  instabilities  in  high-speed  flows  (ASU/DARPA),  and  trajectory  control  of  missiles 
(EglanAFB  /  AFRL) . 

As  a  result  of  this  work,  Eric  Matlis  was  named  a  2002  NASA  Graduate  Fellow.  He  will 
continue  the  work  developed  in  this  grant  in  further  experiments  with  the  two  cone  models 
in  the  Mach  3.5,  and  possibly  a  Mach  6,  “quite”  tunnel  at  NASA  Langley.  Dr.  Meelan 
Choudhari,  in  the  LaRC  Computational  Modeling  and  Simulation  Branch,  will  coordinate  a 
parallel  numerical  study  which  will  be  used  for  experimental  design  and  provide  a  basis  for 
experimental/computational  validation. 


4  Summary 

The  research  program  was  designed  to  add  to  our  understanding  of  transition  to  turbulence 
over  3-D  supersonic/hypersonic  lifting  bodies.  It  brought  together  for  the  first  time,  two 
necessary  capabilities: 

1.  the  ability  to  control  the  steady  and  unsteady  initial  conditions  of  instabilities  leading 
to  transition  to  turbulence,  and 

2.  accurate  time-resolved  measurements  of  flow  quantities  in  compressible  boundary  layers 
using  calibrated  hot-wire  sensors. 

In  addition,  it  teamed  efforts  in  experiments  and  computations,  with  the  computations 
first  guiding  the  design  and  accuracy  of  measurements,  and  experiments  ultimately  providing 
validation  of  the  numerical  simulations. 

Finally,  the  approaches  that  were  developed  for  controlling  initial  conditions  in  our  lab¬ 
oratory  experiments  have  the  potential  to  be  used  in  real  applications  in  transition  control, 
as  well  as  numerous  other  flow  control  applications. 
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